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Abstract. We consider the Mahler measure of the polynomial 1 + x\ + X2 + £3 + £4, which is the first case not yet 
evaluated explicitly. A conjecture due to F. Rodriguez- Villegas represents this Mahler measure as a special value at 
the point 4 of the L-function of a modular modular form of weight 3. We prove that this Mahler measure is equal 
to a linear combination of double L-values of certain meromorphic modular forms of weight 4. 

1. Introduction 



o , 

The logarithmic Mahler measure of a Laurent polynomial 
3' P{x 1 ,...,x n )£C[xf 1 ,...,x± 1 } 
is defined as 

TO ( P ) = fo -\n / log |P(xi,...,a;„)| 

One can show that this integral is always convergent. For a monic polynomial in one variable P € C[x] one can 
E"^ ■ compute m(P) by Jensen's formula 

(1.1) / log|P(z)|- = £ max(0,logH), 

but no explicit formula is known for polynomials in several variables. Let us consider the simplest case of linear 
forms, namely m(l + x\ + • • • + x n ). In 1981 C. Smyth discovered ( Sm ) that 



q /q 

m (1.2) m{l + x 1 +x 2 ) = — L(x- 3 ,2) 

where X _ 3 (n) = (=3), L( X - 3 ,s) = £ = 1 ~ £ + 37 - h + ■ • • and 

n=l n 



7 

1.3) m(l + Xl + x 2 + x 3 ) = r-jC(3). 



, These formulas can be proved by explicit integration. Later we will see another method due to F. Rodriguez- 
£NJ ■ Villegas ( |RV] ) to obtain (|1.2p and (|1.3[) with the help of modular forms. Already in the next case no explicit 
, formula for mil + x\ + x 2 + £3 + £4) is known, and this is the subject of the present paper. One can find in RTV 
■ the numerical value 

m(l + xx + x 2 + x 3 + Xi) = 0.544412561752185... 
and also there is the following conjectural formula. 



Conjecture (F. Rodriguez- Villegas, |BLVD] . see also [Zudlj ): 

'7^15^5 



m(l + xi + x 2 + x 3 + Xi) = 6 ( — — ) i(/i 5 ,4) 

V 2m ) 



where 

fl6 = Viizfvibzf + V (zf7 1 {lhzf = q + q 2 - 3q 3 - 3q 4 
is a CM modular form of weight 3, level 15 and Nebentypus (—7^) ■ 

This modular form arises in |PTV] in the relation to the variety 

{1 + xi + x 2 + x 3 + Xi = 
1 + ^ + ^ + ^ + ^ = 
Xl 12 X3 Xi 

1 
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which can be compactificd to a K3 surface of Picard rank 20. Namely, C.Peters, J. Top, M. van der Vlugt show 
that if X is the minimal resolution of singularities of the above surface then the L-function of H 2 (X) has generic 
Euler factor 

(1 - pT)™ (l - ) P T) 4 (l - A p T + (^) P 2 T 2 ) 

where A p is the pth coefficient in the g-expansion of /15 . 

In order to state our results, consider the modular function t(z) and modular form f(z) of weight 2 

t{z) = _ ( vJ ^) y = - g _ 6g2 _ 21g 3 + „. 



r]{z)r]{3z) 

(V(z) v (3z)r 
{ri(2z)r](6z)y 

for the group 

r (6)+3 = r (6) U {^(^ G SL(2,R) \a,b,c,de z} . 

Throughout the paper we use the differential operator D = = q-^- We need the following modular forms 

of weight 4 

Dt 



.91 = — -/ = l + 2q- Uq 2 + 38g 3 - 142q 4 + 252g 5 - 266g 6 + . . 
t 



(1.5) 92 = — t 91 = ~ q ~ 7q2 ~ 6q3 + 5q4 + U0q5 ' 

£(212i 2 + 251t — 13) , , 

53 = l 9i = 13g + 316g 2 + 2328 9 3 + . . . 

Here g± is indeed a modular form and one can write it as a linear combination of Eisenstein series (see p.8[l ). while 
52 and 53 have poles at the discrete set of points where t(z) — 1. Our main result is the following. 



Theorem. Consider the Chowla-Selberg period for the field K — Q(^f—15) 

(1.6) n 1B = -^(un{- 5 r u Y 4 

3=1 

and the two numbers 

/•OO rOO /"OO /"OO 

(1.7) L(gj,gi,3,l) = (2tt) 4 / gi(is) / / / gj(is 3 ) ds 3 ds 2 dsi ds 

JO Js J si Js2 

for j = 2, 3. One has 

4 

mil + X\ + X2+x 3 + 354) — - mil + xi + x 2 + x 3 ) 

5 



-^-£(flS,ft,3,l) - 1^7^2,51,3,1) 



The reader will find this statement in a slightly different notation under the name of Corollary 16.31 First, let 
us explain why the integrals in (|1.7p converge. For z £ iK+ both t(z) and f(z) are real-valued and one can easily 
check that t(z) < 0. Therefore 52 and 53 have no poles along the imaginary half-axis. When s —> 00 we have 
gi(is) = O(l), 52(25) = 0(e~ 27rs ) and the same holds for 53 because (/-expansions of 52 and 53 start in degree 1, 
therefore the integrated integrals above are convergent at 00. Also one can show that all three functions 9j{is) are 
o(s) when s — > 0, hence they are globally bounded and there is no problem with convergence at s = 0. With the 
help of PARI / CP we find that numerically 

£(52,51, 3,1) = -0.44662442... 

L(53,5i,3,l) = 8.5383217... 

which agrees with the statement of the theorem. 
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These two numbers are periods whose geometric meaning is not clear at the moment. If for example gi were a 
holomorphic cusp form then the number defined in (jl.7l) would be indeed the value of the corresponding double 
L-function L(g2,g%, s%, Sx) at S2 = 3, s\ = 1, which is the motivation for our notation. We discuss double L- 
values of holomorphic modular forms in Section [5] But as soon as forms under consideration have poles in the 
upper half-plane the corresponding multiple integrals become path-dependent and there is no general theory of 
multiple L- values. Also we would like to remark that for the holomorphic modular form g\ in our theorem one has 
m(l + x\ + X2 + £3) = — 5-L(<?i,l), the reader can find the proof of this statement in Section[3J Another observation 
is that the poles of gi ant g% are located at the points from the same field K — (Q)(\/— 15), namely at the images of 
z = g + v/ 24 15 under the group r (6) + 3. 

The structure of the paper is as follows. Sections 2 and 3 follow the approach pioneered by Rodriguez- Villegas 
|RV] . In Section 2 we relate the Mahler measure of 1 + x\ + • • • + x n to the principal period of a pencil of Calabi-Yau 
varieties of dimension n — 1 given by 



and the corresponding Picard-Fuchs differential equation. One needs to do explicitly analytic continuation of 
its solutions from one singular point to another one in order to compute the Mahler measure. For n = 2, 3 the 
differential operators appear to have modular parametrization. This allows us to do necessary analytic continuation 
and derive (|1.2[) and (|1.3[) in Section 3. 

When n = 4 the Picard-Fuchs differential operator is not modular. However one can apply Jensen's formula to 
reduce the number of variables: in Section 4 we observe that in fact m(l + x% + • • • + x n ) can be computed by 
analytic continuation of a solution of a non-homogeneous differential equation with the Picard-Fuchs differential 
operator corresponding to m(l + xi + - ■ - + a; n _i). A non-homogeneous differential equation arises if one considers the 
generating function for the moments of a solution of a homogeneous differential equation along a path. Moreover, 
the differential operator depends only on the initial differential equation being independent of the particular solution 
and the path, while the right-hand side depends on this data (Proposition |4Tj ■ In Section 6 we discuss modular 
interpretation of solutions to a non-homogeneous equation in the case when the differential operator has modular 
parametrization and show that double L-values of modular forms appear naturally in this context. 

Though our main interest is the case n = 4, we keep applying our technique parallelly to the case n = 3 
throughout the paper (Theorem 14.21 and Corollary I6.2j) . This leads to a linear relation (16.41) between a double 
L- value of two Eisenstein series of weight 3, and ordinary L-values L(\3, 2) and C(3). We give a direct proof of this 
relation in Section 7 using a method due to Mathew Rogers and Wadim Zudilin. 

Our original interest in Mahler's measure came from the beautiful paper |RV] which has been inspiring us the 
whole time we were working on this project. We would like to thank our friends and colleagues Sergey Galkin, 
Vasily Golyshev, Anton Mellit, Maxim Smirnov and Wadim Zudilin for their interest in our work. Both authors are 
greatful to the Max-Planck-Institute fur Mathematik in Bonn for providing wonderful working conditions where a 
significant part of this work has been done. 



is well defined for small t since \P\ is bounded on the torus. We call a(t) the principal period of P. It is the 
generating function for the sequence 




2. Mahler Measures and Differential Equations 



For a Laurent polynomial P{x\ 1 . . . 



x n ) the function 




(2.1) 



the constant term of P(xi, . . . , 



•&n) 



m 



since 




00 



2J a m t : 



m=0 
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Suppose that the polynomial P takes only nonnegative real values on the torus {|xj| = 1}. Then the Mahler 
measure m(P) can be computed as follows. For small real t < one has 

= - log(-t) - V — a m = ~(t—) a(t), 
£ — ' m at 

m— 1 

Though we did computation only for small real t < 0, the first integral here and the terminal expression are 
holomorphic in t and defined in some neighbourhood of real negative half-axis (apart from possibly finitely many 
punctures where a(t) has singularities). Therefore 

(2.2) m(P) = -Re (t-£) _1 a(t) 

OX t—oo 

where the analytic continuation is done along — oo < t < and we added real part to be independent of the branch 
of log(t), e.g. we can now assume throughout the paper that 

,00 00 

v at' ' ' tti 

On the other hand, it is known ([SB]) that the sequence (|2.1[) always satisfies a recursion, i.e. a(t) is a solution 
to an ordinary differential equation 

(2.3) C{t,t^ja(t) = 

where £ is a certain polynomial in two non-commuting variables. Finally we see that Mahler measure m(P) can 
be computed by doing analytic continuation of a particular solutions to an ordinary differential equation which one 
constructs from the polynomial P. 

Let us apply this strategy to the linear polynomials. Observe that 

m(l + xi-\ \-x n ) = -m(P n ) 

where 

(2.4) P n = (l + x 1 + ---+x n )(l + — + ••• + — ) 

takes nonnegative real values on the torus. Consider the sequence of the free coefficients of the powers of P n : 

n = 2 a m : 1 , 3 , 15 , 93 , 639 . . . 
n = 3 a m : 1, 4, 28, 256, 2716 ... 
n = 4 a m : 1 , 5 , 45 , 545 , 7885 . . . 
The corresponding differential equations 

£ n (t,t—)a(t) = 

are given by: 

C 2 (t, 0) = 9 2 - t{\Q0 2 + 106» + 3) + 9t 2 (6 + I) 2 

C 3 (t, 6) = 9 3 - 2t(29 + 1)(56> 2 +59 + 2) + QAt 2 {6 + if 



£ 4 (t, 6) = 4 - t(356 4 + 706 3 + 639 2 + 286 + 5) 

+ t 2 (9 + 1) 2 (2596» 2 + 5186* + 285) - 225t 3 (6> + 1) 2 (6> + 2) 2 
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(see |Ver08] for the general form of the operator.) We use the notation 9 = to distinguish it from D = q-£-. 
In all three cases there is a unique analytic at t = solution satisfying a(t) = 1 + o(t). 

The equations ^2(^1,^2) = A and P3,(xi,X2,xgi) = A describe families of elliptic curves and if3-surfaces of rank 
19 respectively. It is therefore natural that the differential equations £2 and £3 have modular parametrization, 
hence we can easily do analytic continuation of their solutions and compute the corresponding Mahler measures by 
formula (I2.2j) . We do this in the next section. 

Unfortunately, this method is not applicable in the case m(l+xi +X2+X3+X4). The equation P^(x\, X2, £3, £4) = 
A describes a family of Calabi-Yau threefolds, and hence we do not expect £4 to have modular parametrization. 
Indeed, one can check that the differential operator £4 is not a symmetric cube of any second order differential 
operator and therefore it does not admit a modular parametrization. Further in the paper we show that one can 
still use the operator £3 to compute m(l + x± + x% + ^3 + £4), though the price paid is that we have to consider 
non-homogeneous differential equations. 

3. Modular parametrizations of £ 2 and £ 3 

Recall that for an arbitrary modular function t{z) and a modular form f(z) of weight k on a congruence subgroup 
of SL(2, Z) one can construct an ordinary differential operator of order + 1 with algebraic coefficients 

fe+i 

(3.1) I>W(!) l < 

i=0 

such that the functions 

(3.2) f(z),zf(z) 7 ...,z k f(z) 
span the kernel of the pull-back of Ctj to the upper half-plane 

X>(*»(^)' - EMh)'- 

t=0 v ' i=Q 

It follows that an operator with these properties is unique up to multiplication by algebraic functions of t on the 
left. On the other hand, the operator 

1 / d \ fe+i 1 
t>(z)-f(z)\d~z) J(zj 
obviously annihilates the local system (|3. 21) and it is a routine to check that if we rewrite it as 

t'(z).f(z)(fc) W) = 

all coefficients gi{z) will be modular functions and hence can be written as some algebraic functions gi{z) = Ci(t(z)). 
The reader could refer to |Zag , Proposition 21 for several constructions of the differential equation satisfied by a 
modular form. 

Recall that D = ^I'h = Ijij' ^ n * ne y i ew °f the above said we make a choice and define the operator 

0.3) ^f = nh Dk+1 r 

This choice corresponds to the leading coefficient in p.l[) being 

D k t 

T 2 

It is not hard to check that both £2 and £3 can be obtained from certain pairs of a modular function and 
modular form, namely the following ones (see [Ver96 for all the details in the case of £3.) 

Proposition 3.1. With 



Cfe+l(<) = , 



(3.4) 



^{Qzf^zf 2 3 

77(3z)V2^ 8 q q q 

ri(2z) 6 ri(3z) „ „ 2 „ , 

f = A3 a d = l + 3g + 3 9 2 + 3g 3 - 
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has 

1 „ / d\ 



Proposition 3.2. With t(z) and f(z) as in (|1.4[) one has 

^ = I* («•*!)• 

Let us use these modular parametrizations to compute m(P2) and m(P^,) by formula (12.21) . With (I3.4[) we see 
that at q = we have t = and / = 1, hence / coincides with a(t) near i = 0. Since t runs over the negative real 
axis when z 6 | + zR= we have by (|2.2D 



m(P 2 ) = -Re ^(t) 



-Re D~ 



T') 



Here and throughout the paper we make a particular choice for D 1 by letting it to act on g-series as 



D 1 = a log q + ^ — 



n 

n— n— 1 



which is in accordance with our choice for (tjj) , so that the above computation is correct. To evaluate the 
terminal expression in (|3.5[) let us recall the definition of L-function of a modular form. 

For a modular form g = X^^Lo c ™9™ °f weight k on a congruence subgroup of SL(2, Z) one has c„ = 0(n k ~ 1 ) and 
the L-function of g is defined by 



H9,s) = 

n=l 

when Res > fc. This function can be continued as a meromorphic function to the whole C. Moreover, if g is a cusp 
form then L(g, s) is holomorphic everywhere in C. 

Proposition 3.3. Let g be a modular form of weight k with Co = and p < k be an arbitrary integer. If L(g, s) 
has no poles with Res > p then one has 

Km (D- p g y q ) = L(g,p). 



If Co ^ the same holds when p < , 



lim g(q) = c + £(fl>,0) 

9->l 



and 



lim (xr^Vg) = lim (c log g + V ^ g») = 1) 



n— 1 

assuming that the branch o/logg is taken so that lim logg = 0. 

The reason we do not consider p > 1 in the latter case is that it is not that clear how to define D~ p when Co ^ 0. 
Note also that one always has L(g,p) = when p < 0. Indeed, the function A(g, s) = r(s)/(27r) s L(<7, s) satisfies a 
functional equation when s goes to k — s and it is obviously holomorphic when Re s > k, hence also when Re s < 0. 
Since T(s) has poles at nonpositive integers L(g, s) has zeros at all integers s < 0. If g is a cusp form then also 
L(g, 0) = by the same reason since A(<?, s) is holomorphic in the entire complex plane. 

Proof of Provosition \3.3[ First let Co = or Co ^ but p < 0. When Re w > k — p one has 

4.W—X „—2nnt 



T(w) . . r(w) x -> C n x ^ C 



n=l 



n 



f" 



ii 



- 1 ( J D-^)(ii)dt. 
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Since T(w) is small when the imaginary part of w is large and L(g , p + w) is uniformly bounded we can apply the 
inverse Mellin transform. Namely, with any real c > k — p one has 

D-*g)(it) = — / -^L(g,p + w)dw 

1 /— e-Moo 

= L(g,p) + — / L{g,p + w)dw 

where we moved the path of integration and < e < 1 is chosen sufficiently small so that C(g, s) still has no poles 
with Re s > p — e. The last integral is 0(t e ) and obviously vanishes when t 0. 

If Co 7^ the calculation above would remain correct if we substitute D~ p (q) by g(q) — cq and D~ l (q) — cq \ogq 
when p = and p — 1 correspondingly. □ 



(3.6) 



Going back to (|3.5[) . consider a modular form of weight 3 given by 

g(z) = (^ff)(z + \) = l + q- r oq 2 + q 3 + llq 4 -24q 5 + ... 
= E 3 , x _ a (z) - 2E 3 , X _ 3 (2z) - 8E 3 . X _ 3 (4z) 
where E 3iX _ 3 € M3(ro(3), X-3) is the Eisenstein series 

(3.7) E 3 , x _ 3 = -i + ££ x _ 3 (d)dV. 

ThenL( J B 3>x _ 3 , S ) = £(«) i(x-3, « - 2) and 

L{9,s) = (l-l-A) C ( S )i(x-3, S -2) 

is holomorphic in the entire complex plane. Since Fourier coefficients of the form g are real L(g, s) takes real values 
at real arguments s. Combining f|3.5[) with Proposition 13.31 we get 

m(P 2 ) - -L(g,l) = 2IimC(«)L(x-3,«-2) = ^L(x_ a ,2). 

Recall that m(l + xi + .t 2 ) = ^m(P 2 ), hence we have just reproved (|1.2[) . 

Analogously, with (jl.4l) we have that t(z) assumes all negative real values along the imaginary half-axis and 
t = 00 at z — 0, hence 



m(P 3 ) = -Re D- 1 (^/ 



9=1 



We consider 



g = —f = l + 2q- Uq 2 + 38q 3 - U2q A + 2h2q 5 - 266g 6 + . . . 
(3.8) t 

= 2E 4 (z) - 32E 4 (2z) - 18E 4 (3z) + 288P 4 (6z) 

with the Eisenstein series 

n>l d\n 

The function L(E 4 , s) — ((s)C t (s — 3) has the only pole at s = 4, and one can easily see that 

/ x / 32 18 288\ 

L{g,s) = (2---- + — )c(*)C(«-3) 

is holomorphic in the entire complex plane because the factor in the brackets vanishes at s — 4. Finally, 

- -x (9 ,D - -^ + ^)ii s cWC(»-3) - ™. 

This again reproves (11.31) because m(l + xi + .X2 +^3) = ^^(Ps). 
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4. Computation of m(P n ) via £„_i 

Observe that in general the Mahler measure of a Laurent polynomial P which takes nonnegative real values on 
the torus {\x\\ = ■ ■ ■ = \x n \ = 1} can be written as 

m{P) = J miX log(A)a*(A)dA 

where 

A min = min P(xi, . . . ,x n ) , 
|a;i|=-=|a„|=l 

A max = max P(xi ,...,x n ), 

\xi\=-=\x n \=l 

and a* (A) is equal to the integral over the variety 

{P(xi,...,x n ) = A}n{|x x | =••• = |x„| = 1} 

of the the (n — l)-form lu\ defined as the residue 

1 dx\ dx n 

A ... a = u\ A dX , 



{2iti) n X\ x 

i.e. 



a* (A) = / uo x . 
J{p=\}n{\x t \=i} 

Along the same lines as we did in Section^ one can recover m(P) form the generating function of the "moments" 
of a* (A) 

a m = / A m a*(A)dA, a(t) = V a m t m 



m=0 



by formula (|2.2p . Indeed, the moment a m is exactly the constant term of P m , so this a(t) is identical with the one 
in the previous section. Also one can see it directly by repeating the old trick in our new notation: for — oo < t < 

log(A- -) a*(X)d\ 

00 J- 771 /*A max 

= -log(-i)- — / X m a*(X)dX 

m=l "' A min 

= - log(-t) - V — a m = a{t), 
£ — ' m at 

m—l 

hence 

r-Amax 



m(P) = I \og(X)a* (X)dX = -Re (t^-) 1 a{t) 

JXmin tit 

where we again assume that the analytic continuation of ait) is done along the real negative halfaxis. In Section[5]we 
mentioned without a proof that a{t) satisfies a differential equation (|2.3p . Now this fact follows from the proposition 
below because a* (A) is a period for the 1-parametric family of varieties P(xi, . . . ,x n ) = X and it satisfies a certain 
differential equation 

£(A,A^)a*(A) = 0, 

namely the Picard-Fuchs differential equation for this family. Proposition 14. 1 1 states that moments of a solution of 
a differential equation satisfy another differential equation determined by the initial one, though in general they are 
solutions of this equation with a right-hand side. This right-hand side is a simple rational function which depends 
on the path of integration and on the choice of the solution along this path. In the above situation with a* (A) 
along the real line from A m i n to A max the right-hand side exceptionally appears to vanish. But later we will need a 
general case as well. 
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For a Laurent series J2 n a ^ n we introduce the notations 

n 

[E a «* 



n>0 
n<0 



for its parts with nonnegative and negative powers correspondingly. One then has 

J2 a nt n = [E a « fl L + [E a « fl ] + - 

n n n 

Proposition 4.1. For a polynomial differential operator 



M N , 

C(\6) = EE C ^< 6 = X dX 

i=0 i=0 



consider a solution F(X) of CF — along some path between X — a and A = j3. Then for the generating function 
of its moments 



b(t) = E r / xnp ( x ) dx 



one has 



c(±,-6 t -l)b(t) = h(t) 

where the right-hand h(t) is a rational function which can have poles at most att = 0, ^, j and is defined as follows. 
Let 

M N 

c^{\,e) = EE c « A< ^" fe 

i=0 j=k 

and for given A consider a rational function of t 



N-l 



H x (t) = A E(^'F)(A)[£« +1 )(± -Ot-l)^ 

3=0 



Xt 



The 



h(t) = -0 t -l)b(tj\ - H fj (t) + H a (t). 



Proof. Integration by parts yields 



/ \ i (e : >F){X)d\ = Y, x (-i- 1 )* F ( x ) + (-i-iy \ i F{\)d\ 

Ja s=0 *=« J a 
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and we apply this formula to every term below to get 



ij n—0 



= E*" / X l U(^0x)F)Wd\ = E^E*" / A n+i (0 J >)(A)dA 

n=0 " — " 

oo i-i 

E' n E%-[E^ +1 ( 



A=/3 



n—0 i.j s—0 
N-l 



-7^-^-l) s 6l J - 1 - s F(A) + (-n - i - l) 3 / A n+, F(A)dA 

A— a 



^^•f(a) ^E A " + **" +i ( 



fc=0 



j>fc+l n=0 



iy-i-h 



A=a 



1 °° /-/3 

E c ^ 77 E - 1 l ) Hn+i / A" +l ^(A)dA 

i,j n=0 ,/a 



- H a (t) + [£(-, -0 t -l)&(t) 



) 1 




A=/3 

+ 


£(i -0 t -l)&(t)' 


/1-At 


+ 


A=a 


-l)&(t)' 


+ 







□ 



Now we introduce another idea which will allow to apply the above proposition to our case of linear polynomials. 
We consider P = P n as defined in (|2.4p . and let C n be the corresponding Picard-Fuchs differential operator. It can 

be easily recovered from C n since (up to a simple multiplier) the operators C n (t,9t) and C n (h, —6 t — lj must be 
equal. For example, with 

C 2 = 96» 2 -A(lO6» 2 + lO6i + 3) + A 2 (0 + l) 2 
C 3 = 64<9 3 -2A(2(9 + l)(56i 2 + ! 



2) + A 2 (6> + l) 3 



1 

: 9^ 
1 



C 3 (64t,6] 



t 2 



one can easily check that 

(4.1) Z 2 (i-< 

and 

(4.2) Hi'- ' 
Let a* (A) be the solution of C n a* = such that 

f(n+l) 2 

(4.3) m(P n ) = / log(A)a*(A)dA. 

Jo 

Applying Jensen's formula (jl.ip in the variable x n +i gives 

2 m(P " +l) = (2^T J / N|=1 l0g|1 + ' Tl + ''- + a; " +l1 — 



dx\ dx n +i 



(2tt»)» 4 



log 1 1 + El 



or 



(4.4) 



|=l,|l+iBiH hxnl>l 

. (n+lr 

m(P nH 



^n+l 
X\ Xji 



log(A)a*(A)dA. 



Observe that (14.31) and (|4.4I) differ only by the lower limit of integration. This approach yields us to state that for 
every n there is a simple rational function h n (t) and an analytic solution b n (t) of 

£ n (t,Ab n {t) = h n {t) 



such that 



m(P n+1 ) = -Re (t— ) 1 6„(t) 
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Below we give exact statements for n — 2 and n — 3. We will formulate our results using solutions with rational 
coefficients rather then their transcendental linear combinations. This makes sense from the number-theoretical 
point of view and will be used later. 

Theorem 4.2. Take C 2 (t, 0) = 6 2 - t(109 2 + 106 + 3) + 9t 2 (0 + l) 2 and consider the following analytic at t = 
solutions of 

C 2 (t,t—) 4>{t) = (j){t) = l + 3t + ... 

c 2 (t,t-)m = -^— t m = t+... 

Then 

d \-i r 3 6 
~dt> 



m(P 3 ) = -Re(i-) [-<j>(t) + —ip(t)] 



Though all Mahler measures in this theorem were already computed before, we will use our result to relate them 
to double L- values of modular forms later. Next theorem already deals with the interesting case n = 4. 



Theorem 4.3. Let fiis be the Chowla-Selberg period for the field K = Q(v / — 15) as in (|1.6|) and b{t) be the unique 
analytic at t = solution of the non-homogeneous differential equation 

rU+ d \uM 3^^ i(212i 2 + 25K - 13) 3^5 t 
C 3 (t,t-)b(t) — + 

satisfying b(t) = 5 + 0(t). Then 

m(P 4 ) = -Re^)- 1 ^) 
5. Proofs 

Recall from Section [5] that the function 

If 1 dx\ dx r< 



a ^ (27ri)™ J\ Xl \^... = \ Xn \^ l 1 -tP n (xi,...,x n ) xi 
is the unique analytic at t = solution of C n a = satisfying a(0) = 1. Let us write a(t) as 

f in+lf «*W JX 



1-tX 

where a* (A) is a solution of C n a* = 0. Let us also introduce 

b(t) = I r^x dX = bo + blt 

f 1 o*(A) „ 
c(t) = J ^—jk = °o + cit + ... 

In Section @] we showed that 

m(P n+1 ) = -Re 

In this section for n = 2, 3 we will compute the coefficient bo and prove that 

h(t) := £„&(*) = -C n c(t) 

has the form as in Theorems 14.21 and 14.31 correspondingly. The proofs of Theorems 14.21 and 14.31 go along the same 
lines. First, we identify the solution a* (A) for A G [0, 1] in terms of the Frobenius basis in the space of solutions 
near A = 0. In order to do this we use asymptotics of a(t) when t is large, and we use the modular parametrization 
of the differential equation to find this asymptotics. As soon as we know a* (A) explicitly, we apply Proposition 14. II 
to finish the proof. 

We note that the differential operator £2 has singularities at A = 0, 1 which are regular singular points of maximal 
unipotent monodromy, whereas the operator £3 has a regular singular point of maximal unipotent monodromy at 



12 EVGENY SHINDER, MASHA VLASENKO 



A = and is nonsingular at A = 1. The period Q15 appears in Theorem 14.31 because the point A = 1 corresponds 
under the modular parametrization of £3 to a CM point of conductor 15. 
It hopefully will not confuse the reader that we use the same notation 

a(t),a*(X),b(t),c(t) 

for the integrals corresponding to the case n = 2 in Lemmas 5.1 — 5.3 and to the case n = 3 in Lemmas 5.4 — 5.7. 
Proof of Theorem \4-S\ Here a(t) = </>(t). From Proposition 14. II and Lemmas 15.21 and 15.31 below it follows that 



A ' / 7T 2 VI - t tJ 1T 2 t(l-t)' 

According to (|4.ip we then have 

6 t 



C 2 (t,9)c(t) = t 2 £ 2 (j, -0-l)c(t) 



7T 2 1 -t 



and then 

c 2 (t,e)b(t) = -c 2 (t,e)c(t) (l ' 



7T 2 1 -t 

From Lemma 15731 we have 60 = 1 — Co = f and therefore b(t) = j<fi(t) + -^ifj(t). To finish the proof it remains 
to verify the three lemmas that follow below. 

□ 

Lemma 5.1. 

a* (A) = -4- + 0(A) 

V07T A^O 

a*(A) = -J-log(l-A) + 0(l) 
471-2 v+i- 

Proof. Using the modular parametrization (j3T~4j) we find that when t — > — 00 along negative real axis (this corre- 
sponds to z going down the ray ~ + zR + ) 

(5.1) ta(t) = -L-log(-i) + 0(1) 
On the other hand 

/ \ , , f 9 a* (A) „ f 9 a* (A) „ 

(5.2) ta(t) = / v \ riA = - / — -fdA 



/ 1/t — A J s + A 

Let us write a* (A) = ag^o(A) + ax0i(A) in terms of the Frobenius basis 

0o(A) = 1 + 0(A) 
01 (A) = log(A)0 o (A) + O(A) 
of solutions near A = 0. One can easily check that for any e > 

dX 



o s + A 



logs + 0(1) 



logAdA 1 n _ 2 



S + A -2( l0gS ' + °W 

and comparing (|5.1[) and (|5.2|) we see that ot\ = and ao = ^j^. 
Now let 

«o(A) = 1 + <3(A-1) 

ki(A) = log(A - 1)k (A) + 0(A - 1) 

be the Frobenius basis at A = 1 and a* (A) = ao«o(A) + ai/ii(A) when A — > 1_. Using our modular parametrization 
we find that 

a* (A) 1 f(z) 3 

Ofi = urn ; = urn — = ; — - = . 

1 " ' ...... ' s 47T 2 

□ 



A-n_log(l - A) s^o+ ^ log(l - 9t(z)) 
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Lemma 5.2. In the notation of Proposition ^ ■ 1\ applied to C = £2, a = 0, j3 = 1, F(X) — a*(\) we have: 

H (t) = 

Proof. In the course of the proof we rely on the asymptotics given in Lemma l5.ll Since A = 0, 1 are singular points 
we are going to compute Hq (f ) and H\ (t) as the corresponding limits of H\ (t) . We have 



£ 



(2) 



a- 



and 



£ 



1 U' 



1-Xt 



1 



1 



n 10 1 

q 1 

t t 2 J 1 - Xt 



9t 2 - lOt + 1 1 A - 10 



(A-l)(A-9) 



t 2 (l - Xt) 



t 2 



t 



1 - Xt 



-(9 



1 - Xt 
10 1 



(9-y + ^)M-D + (-y + 2t 2 ; 



1 - xt 



10 



2t 2 



1 



1 



(A-l)(A-9) 



t t 2 '(l-Xt) 2 y t 'l-Xt t 2 (1-At) 2 

These functions have finite limits when A — > and now we see that Ho(t) = because a* (A) is analytic at A = 
and therefore lim A 6 3 'a* (A) = for any j > 0. 

A->0 



Since lim (A - 1) a* (A) = and lim (A - 1) 0a* (A) 



A-s-1 



A^l_ 



we find that Hi (t) = ^rf 



(l-t) ■ 



□ 



Lemma 5.3. The first coefficients cq in the power series expansion of c(t) is equal to \ and we have 



-6- l)c(t) 



Proof. It is easy to compute that 



6 

~^~t 
~3c + ci 



Using the modular parametrization (|3.4p (with modular * and / from (|3.4[) and A = 9t one has £\j = gj£2(A, A ^3-)) 
we compute that 



and 



Co 



a*(A)dA 



9 



v^tt 
18 



f(z)t'(z)dz 



jj=j f{zft{z){l~9t{z))(l-t{z)) 



ds = — 

z—is 4 



c\ — 3co = / (A — 3)a*(X)dX 



18 



^= y (9t(z) - 3)f(z)h(z)(l - 9t(*))(l - t{z)) 
6 



r/.s 



□ 



Proof of Theorem \4.3\ We let Q — ^15. According to Proposition 14. II and Lemmas 15.51 and 15.61 below we compute 
that 

3a/50 2 /13 -13i 2 + 251t + 212\ 3^5/1 1 



-e-x)c(t) 



10 



7rAT 



(i-*) a 



3V5fl 2 — 212i 2 - 251* + 13 3%/5 1 



10tt 



t(l-i) 3 



5tt 3 ^ 2 t(l-t) 



14 EVGENY SHINDER, MASHA VLASENKO 

Now it follows from (14.21) that we have 



-h(t)=L s (t,6)c(t) = -i 2 £ 3 (^-#-l)c(i) 

Sy/Stt 2 t(212t 2 + 251t- 13) 3-y/B" i 



Therefore the function 

" ,6 a* (A) 



satisfies £36 = hit) and its power series expansion at t = starts with 6o = ao — Co = |. The proof will be 
finished after verifying the three lemmas below. □ 

Lemma 5.4. In terms of the Frobenius basis 

<MA) = 1 + 0(A) 

0x(A) - log(A)0 o (A)+O(A) = log(A)+ (l) 
02(A) = log(A) 2 o (A) + O(A) = log(A) 2 + o(l) 

of solutions near A = we have 



At A = 1 



a *( A ) = -g^faiW ~ 61o § 2 ^(A)). 
o*(l) = 0.1649669005300320... = ft 2 (fi = fii 5 ) 

27T 

3 \J% 

6»a*(l) = -0.032993380106006... = — Q 2 

107T 

6>V(1) = 0.00330836512971504... = ^(H^ L). 

150 V 7T 7T J i2^/ 

Proof. With the help of modular parametrization (jl.4p we find that when t — > — oo along negative real axis (this 
corresponds to z going down to along the imaginary axis) 

Mt) + I ^iog(-4) 2 -o 

On the other hand 

" 16 a *( A ) ^ _ f 16 a*i\) 



ta ® = I W^x dX = -J ~s + \ dX 



and since for any e > one has when s —> 

dX 



o s + X 



-logs + 0(1) 



logAdA 1 o 

- -2(log S ) 2 + 0(l) 

^f^-4(log S ) 3 + 0(log (s ) 2 ) 

3 9 
we find that c*2 = 0, ai = — 7-^7 and ao = - — 7 log 2. 

87T^ Att z 

We indicate how to find the values J 'a*(l) for j — 0, 1, 2. Again we use modular parametrization (|1.4j) . namely 
with modular i and / from (11.41) and A = 64i one has = 64T^ 3 ( A ' A 3x)- This \(z) takes real values from the 
interval (0, 1] for z — \ + is and s G (+00, (let us denote r = | + ^ g 15 so that A(r) = 1). Using asymptotics 
at 00 one can check that on the vertical half-line from r to 00 



Now the problem is reduced to computing the values of modular forms and their derivatives at a CM-point of con- 
ductor 15, this leading to expressions involving 51 and tt, see |Zag| , Propositions 26, 27 and Corollary of Proposition 
27. 
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Lemma 5.5. In the notation of Proposition applied to C — £3, a — 0, j3 — 1, F{\) = a*(X) we have: 

ff (t) = 

3fi 2 v / 5-13t 2 + 251i + 212 3^/5 1 



Hi{t) 



IOtt 



(l-*) £ 



5vr 3 ft 2 1 - t 



Proof. One easily checks that Ymv\6 3 a* (A) = 0, whence Ho(t) = 0. In order to compute Hi(t) we need 

A— >o 



C 



(3) 



(f 



-Ot-1 
-Ot-l 

-Ot-1 



1 - 




45 


1 -t- 


+ 


1 -t 


1 - 




9(t-6) 


1 -t- 


+ 




1 - 




29 + 68< - 7t 2 


1 -t- 


+ 


(1-i) 3 



and then the formula for f/i (i) follows after a simple computation with values of O^a* (1) provided by Lemma I5~4l 
Lemma 5.6. TTie /irsi coefficient cq in the power series expansion of c(t) is equal to ^ and we have 

1 „ .\ ,.J /39V5„2 , 3^/5 \ 1 



Proof. It is easy to compute that 



£ 3 (-,-0-l) c (t) 



-ft 2 

V IOtt ' bir^ 2 J f 
4c - ci 



3 f T 1 

<z*(A)dA = — ——, r ■ 2iri / (z — —) f(z) X'(z)dz 
8tH 7 loo v 2 y 

here we use that A = 64i and (g— ) / f 2 = < 2 (1 - 4i)(l - 16*) 



■ (2™) 2 • 64 / (z - -) f(z) 2 t(z)y/(l-4t(z))(l-mt(z))dz 



= 96 / a f(z) 2 t(z) y/(l - 4i(z))(l - 16t(z)) 

' VIE 



ds = — 

z — ^+is 5 



and 



Cl -4c = / (A-4)a*(A)dA 



96/ s/(z) 2 (64i(z) - A)t(z)y/(1 - 4i(z))(l - I6t(z)) 



= -0.708951451918989714. 
^ 2 



7a*(l) -9 6>a*(l) +45 6>V(1) 



ds 



39-\/5 2 Sv 7 ^ 



IOtt 



5tt 3 ^ 2 



6. Double L-values of modular forms 
In Theorem 14.21 we are led to evaluation of 



where ip is the unique analytic at t = solution of 



On . . 
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which satisfies the condition ip(t) =t + o(t). The same happens in Theorem 14.31 Putting this situation into a more 
general context, consider a solution of a non- homogeneous differential equation of order k + 1 which has a modular 
parametrization, i.e. 

C tJ $ = h(t) , 

where t is a modular function, / is a modular form of weight fc, £ t .f is defined by (|3.3p and h(t) is a function of t 
which will be just rational in our cases. We can consider ip — ip(t(z)) as a function in the upper half-plane and we 
rewrite the above differential equation as 

1 :^ fe+1 T = h (t), 



or 



Dt-f f 



£ fc+i y = h{t)-Dt-f. 

Therefore — is an Eichler integral of the modular form h(t) -Dt-f of weight k + 2. (This conclusion is precisely 

the statement of Lemma 1 in |Yanj .) Let us assume in addition that the modular function t takes values and oo 
at q = and q = 1 correspondingly. Then 



( t l)-Vw| (TO = D -(^)| a=i = D -'( 92D -'-v) 

where 

Si = h(t)-Dt-f, 92 = 

are two modular forms of weight fc + 2. According to the proposition below the above expression appears to be the 
double L-value of these two forms. Let us give the definition. 

Let gi — J2 n >o a nQ n an d 92 = J2 m >o bmO m be two modular forms of weight fc on a congruence subgroup of 
SL(2, Z), and let in addition a = 0. Their double L-function (denoted by L* in [5r] ) is defined for Re(si + S2) > 2fc, 
Res2 > k by 



Dt 



EV" 1 On^m 
n Sl (n + m) s 

n=lm=0 y ' 



The question of the analytic continuation simultaneously in the two variables s±, S2 is rather tricky and we do not 
want to consider it here. But it appears that for any fixed integer s\ = p the function L(gi, 32, S2) is well-defined 
for S2 with sufficiently large real part, and we can easily prove analytic continuation in this variable when p > 0. 
In order to do this one writes ([5r]) 

p-i 



(6-1) L(g 1 ,g 2 ,p,s 2 ) = r ^ r ( g2 ) H9u92,P - m, s 2 + m) 

where 



/•OO />CJO 

A(g 1 ,g 2 ,s 1 ,s 2 ) = J t S2 - 1 g 2 (it) J v Sl ~ 1 g 1 (iv)dv dt . 
Now observe that these integrals are well defined for all Si,s 2 . Indeed, it follows from the estimates 

/oo 
« ai_1 5i(*t;)d« = 0(t si - 1 e- 2wt ) , i->oo 

= 0(t si - k e-^), t^Q 

since gi(it) = o(t~ k e~ 2 ?~^J when < — > 0. Therefore formula (|6.1[) gives analytic continuation of L(gi,g 2 ,p, s 2 ) in the 

variable S2 with integer p > 1. Moreover, this function is holomorphic in the entire complex plane because 1/T(s 2 ) 
is holomorphic, and we can speak of "double L-values" L(gi, g 2 ,pi,p 2 ) with integers pi,P2 whenever p\ > 0. Notice 
also that L(g±, g 2 ,Pi,p 2 ) = if p 2 < as one can see from (|6.ip since T(s2) has poles at nonpositive integers. 

Proposition 6.1. Let gi,g 2 be two modular forms of weight k on a congruence subgroup, g\ vanishing at 00. Then 
for any integers < p\ < k and p 2 < k one has 



(6.2) UmD-»(fla.D-"^i)(g) = L( gi ,g 2 , Pl ,p 2 ) 
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Proof. Whenever Re (p 2 + w) > k one has 

T(w) T(w) ^ ^ a n b m 

K ' v ' n=lm=0 V ; 



CO OO , />00 

^ nPi (n + m)P* J n 

n=lm=0 V /JO 

J™ t w ~ l D~^ (g 2 • I?- pi . 9 i) (it)dt . 



By Mellin's inversion theorem with an arbitrary real c > k — P2 one has 

D^( 92 .D^ 9l )(it) = ^-j c ioo J^L( 91 , 92 , P1 , P2+W)dw 

1 /■-*+<«> T r w \ 
= L(gi,92,Pi,P2) + ^— J w L{g u g 2 ,p\,P2 + w)dw 

with any < e < 1 and we moved the path using the fact that L{ 9 \, g 2 ,p\ 1 p 2 + w) is everywhere holomorphic in 
w. The last integral obviously vanishes when t — > 0, and (|6.2[) follows. □ 

In the case of Theorem 14.21 we use the modular parametrization (13.4[) and h(t) = j^-r. Since t — oo at z = ^ we 
consider the shifted forms 

9l (z) = = l + q-5q 2 + q 3 + llq 4 -24q 5 + ... 

= E 3 , x _ 3 (z) - 2E 3tX _ 3 (2z) - 8£ 3 , x _ 3 (42) 

(6 - 3) / £>i \ 1 

52(2) = (7— ^ = -9 - V - 9 3 + 16 9 4 + 24g 5 - 4 9 6 + . . . 

= -£ 3 , x _3 (z) - 7E 3 . X _ 3 (2z) + 8£ 3 , x _ 3 (4^) 

where E 3tX is the Eisenstein series defined in (|3.7p . The form 51 already appeared in (|3.6D and we had that 
m(P 2 ) = —L( 9 i, 1). Using Proposition 16 . 1 1 we now rewrite the statement of Theorem 14.21 as follows. 

Corollary 6.2. With the modular forms g\,g 2 of weight 3 defined in (I6.3P one has 

m(P 3 ) - -MP 2 ) = -^L( 92 , 9l ,2,l). 

4 7T^ 

Plugging in the values of m(P 2 ) and m(P 3 ) which we compute from (|1.2[) . f) 1 . 3[) into the formula given in Corollary 
we obtain the following relation between double and ordinary L- values: 



(6.4) L(g 2 , g lt 2, 1) = ^L( X _ 3 , 2) - |c(3). 

We give a straightforward proof of this relation in the next section. 

For the Theorem 14.31 we use the modular parametrization (ll.4j) and we have to consider two solutions with 
h(t) = and h(t) = 212t +^ 3 t ~ 13 _ Also we have t = at z = ioo and t — 00 at z — 0. According to our strategy, 
we define the modular forms of weight 4 

9l = —f = l + 2q- Uq 2 + 38q 3 - 142q 4 + 252q 5 - 266g 6 + . . . 

.92 = YZrJ = -9 - V - 6q 3 + 5g 4 + 120g 5 + 498g 6 + . . . 
212t 2 + 251t- 13 

(observe that they are the same ones as in (|1.5p ). Here 171 is a holomorphic modular form, it already appeared 
in (13.8[) and there we had m(P 3 ) = —L(gi, 1). The forms g 2 and g 3 are meromorphic with the poles at the points 
where t = 1. Using the fact that t has no poles on imaginary half-axis we defined the corresponding double L-values 
L(g 2 ,gi,3,l), £(53,31, 3,1) in (TE7]) in Section[U 



•93 = T, Ixa Dt 'f = 13<7 + 316 9 2 + 2328g 3 
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Corollary 6.3. With the double L-values defined in (11.71) one has 



^1£ 

We know from Proposition I3.2I that this 



m(P 4 ) - -m(P 3 ) - v 107r 1J L(g 3 ,gi,3,l) - ^^r L ^gi, 3, 1) 



t— oo 



Froo/. Due to Theorem [O] we have that m(P 4 ) = -Re (i^) 1 6(t) 
differential equation has modular parametrization by t(z) and /(z). Therefore 

*«•» - f/M + /Wi>-(-^»W + i^-w). 

We use the path in the upper halfplane from z — ioo to z = along imaginary half-axis, exactly where one has 
— oo < t(z) < 0. As it was explained in Section [TJ all three forms are holomorphic along this path. For gj(z) with 
both j = 2, 3 we then have 



D -i, 



pOO pOO pOO pOO 

[gi- D~ 3 gj)(iv) = (2n) i J gi(is) J J J gj(is3)ds 3 ds 2 ds 1 ds 

and therefore the numbers (|1.7|l are the limiting values at v = 0. □ 

7. Explicit computation of a double L-value in formula 
In this section we show how to compute the iterated integral 



(7.1) 



r°° i r° da 

/ 9l (z)-D- 2 g 2 (z)dz=— gi ( q ).D~ 2 g 2 (q)^ 
Jo 2tti Jx q 



-^-.D-\ 9l -D- 2 g 2 ) 
2m 



9=1 

for two modular forms g\ ,g 2 of weight 3 defined in (|6.3p which leads to an alternative proof of formula (|6.4D . 

We use a powerful method due to Mathew Rogers and Wadim Zudilin [RZT] . [RZ2] . [Zud2] . [Zud3] of computing 
double i-values of Eisenstein-like series. We are grateful to Wadim for explaining to us his method and it's 
applicability in this situation. Unfortunately, the more complicated and exciting L-values from Corrollary 16.31 do 
not seem to be computable in the same way due to lack of Eisenstein-like representation for the forms (71,52,53 
refered to in Corollary 16.31 

We briefly describe the method as follows: the Atkin-Lehner involution z — > — is applied to gi, the resulting 
modular form being denoted by gi(z): 

gi(z) = const ■ 9i(-^) z ~ 3 , 
so that the integral (|7.1I) will take the following form: 



const - / g~i{—T7r)D g 2 {z)z dz. 
Jo L2z 

We then expand the integral as a quadruple sum, make a variable change and collapse the sum back in order to 
get 

/>ioo -t 

const - / fi(u)( const + f 2 ( — ))udu, 

Jo 12u 
where f\ and f 2 are Eisenstein series of weight 1. We apply the Atkin-Lehner involution again: 

h{-T^-) = const ■ h{u)u, 

12u 

this time rewriting the integral as 

pioo 

const - / fi(u)(const + f 2 (u) u)u du — const ■ L(fi,2) + const ■ L(fif 2 ,3), 
Jo 

with fx an Eisenstein series of weight 1 and character x-3 yielding £(2) • L(xs, 2) and f\ f 2 an Eisenstein series for 
SL 2 {Z) of weight 2 yielding C(2) • C(3). 

Note that no regularization is neccessary in our integral, since g 2 vanishes at z = 00 and <?i vanishes at z = 0. 
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In the course of the computation we will need to apply Atkin-Lehner involution to Eisenstein series. For that 
we express Eisenstein series as linear combinations of eta-products. Let N > 1. Consider an eta-product 



f(z) = Y[v(d 



3 

3 



where kj are integers and dj non- negative integers dividing N. Let d'j = N/dj and 

f(z) = Y[r 1 (di-z) k >. 

j 

We have 

(7.2) /(-i-) = {^rz w \{df /2 f{z) 

3 

for w = \ . kj (the weight). This follows from the basic transformation formula: 



i]( — ) = V—izrj(z). 



We use the following two Eisenstein series of weight 3 

x = -n^-TT =~l+ E X-3(n)n 2 q nm = -^ + q -3q 2 + q 3 

i>i 

E 3 , X -s^)^ !I ^T= E X~3(rn)n 2 q nm =q + 3q 2 + 9q 3 + 13q 4 + ... 



We apply (j7T2")) : 



1 V {zf 1 
9?7(3z) 3 9 


r?(3z) 9 
r](z) 3 


E X- 

n , m > 1 


9i{z) = [1 


,-2,-8]- 


92{z) = [- 


-1,-7,8]- 


#3, X _ 3 0) 


= £ 

i 

12z 




= E 

i 

12s 




12z 



:X- 3 ( z ), £; 3, X -3( 2z ),^3, X - 



• X - 3V 12z 
• x - 3(_ 3l 



_\ _ _o3o5/2- 3 j 



and hence 

[a,b,c] • [^.sW.VsfS^VsMK— ) = -3 5 / 2 i2 3 [ C! 2 3 6,2 6 a] • [E 3 , x _ 3 (z),E 3 , x _ 3 (2z) , E 3 , x _ 3 (4z)}. 
In particular for gi(z) = [1,-2,-8] • [E 3-x _ 3 (z), E 3:X _ 3 (2z), E 3 _ x _ 3 (4z)] we have 

.9i(" T ^) = 8-3 5 / 2 zz 3 g 1 (z) 

= [1,2,-8] ■ [£ 3 , x _ 3 ,£«_ 3 ,£g_J(z), 
= q + 5q 2 + 9q 3 + llq 4 + 24q 5 + ... 

or equivalently 

iz~ 3 1 

(7-3) = -_ 7 _ §l( __) 

Formula (I7.3[) allows us to rewrite the iterated integral (|7.1[) as 
(7.4) -VSmL h{-^)E- 2 ^z)z~ 3 dz 



20 



EVGENY SHINDER, MASHA VLASENKO 



We now make use of quadruple sums. For that we write our form g\ and g 2 as: 

gi{z)= J2 ai(m 1 )b 1 {n 1 )nlq mini = q + 5q 2 + 9q 3 + llq 4 + 24q 5 + . 

mi ,ni> 1 

g 2 (z) = a 2 (m 2 )b 2 (n 2 )n 2 2 q m ^ = -q - Aq 2 - q 3 + l6q A + 24q 5 + 



m 2 ,n 2 >l 



oi(m) = x-3(m) 

bi(n) = 1 + — [n even] — — [n divisible by 4] 
a 2 (m) = — 1 — 7[m even] + 8 [to divisible by 4] 
h(n) = X-a(n) 
We obtain a quadruple sum by plugging in the expansions 

mi ,ni>l 

D~ 2 g 2 (z) = ^ a 2 (m 2 )b 2 (n 2 )^exp{2^im 2 n 2 z) 



TOS 



into flLH) : 

1 



233V2 7 §i(- I § i )^- 2 ff 2 (^)^ 3 ^ 

2331/2 ^ ai(mi)6i(ni)a 2 (m 2 )fe2(n2)^ y exp^27ri( — j^— + m 2 n 2 z)jz~ 3 dz. 



1 JQ 

mi ,m ,7712 ,ri2 



Then we change variables. First we let w = — and obtain: 

1 o2 ■ 2 pioo 

12 % ^-^ nf / / m 2 n 2y 



, a 1 (mi)b 1 (ni)a 2 (m2)b 2 {n2)^ expy2iri{miniw fj^f)) 7 "^ 

mi, Til ,7712 ; n 2 ^ 



and then w = : 

m 2 



2 • 3 3 ^ 2 i di(TOi)6i(ni)a2(TO2)&2(^2) / ex-pi 2 r Ki(m\vn 2 u — -)judu 



mi ,ni ,m2 ,n2 

which we collapse to 

/■ io ° 1 1 

(7.5) ~ 2 ' 33/ 7 ^( u )(^(-i^)-6)« du 

where 

h(z)= ]T ^i(mi)a 2 (m 2 )q m ^ 

mi ,m2>l 



6 

ni ,n2> 1 



(the term | makes f 2 (z) a modular form. 
Let 



E 

m,n>l 

then 



m,n>l 

%_ 3 (m)[n divisible by l]q nm = E\{1 z) — -, 
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so that 

fa(z) = -Et(z) - 7E 1 (2z) + &Ei(4«) = -q - 7q 2 - q 3 + 7q 4 - 7q 6 - 2q 7 + . . . 
fa(z) = E x {z) + ^E 1 (2z) - \E^z) = 1 + q + l/2q 2 + q 3 + l/2q 4 + l/2q 6 + 2q 7 + 



In order to evaluate (|7.5|1 we apply Atkin-Lehner involution to fa. We first express fa as a combination of 
eta-products: 

1 f](4z) 2 r](12z) 2 1 i](2z) 6 T](3z) 



and then use 



^ 2 ^ 2 77(22)77(6,2) ' 6 i](z) 3 r)(6z) 2 
l ( 32 V /2 / ■ sViZzfriiz) 2 If 6 6 4 \i/ 2 . v (6z) 6 r](4:z) 



f( _J_ ] = l(J^\ 1/2 (_- M3zYy{zy 1/ 6 b 4 nV2 
721 12z j 2v6-2/ 1 * Z \(6z)r ? (2z) 6V12 3 2V 1 ^ Vl2z) 3 ?7(2z) 2 



, 1 7?(3^) 2 77(63)677(4*) _ . 

/2(Z) " 2 7,(6z)77(2z) + ?7 (12,)3^(2z)2 " + 2i?l(2z) + 8i?l(4z) 



= 3/2 - g + 2g 2 - q 3 + 7q 4 + 2q 6 + 
We continue rewriting the integral (|7.5p using 



(fc-1)! 



!{z)z dz = F2# L(/ ' fc) 



3 1/2 i / /i(u)udu-3 2 / fa{u)fa{u)u 2 du 



3 1 / 2 ?' 3 2 • 2 

r L(/ 1 ,2)-- r —^L(/ 1 (u)/ 3 («),3) 



-27Ti) 2 W " ' (-27Ti) 3 

Z 1/2 i , 3 2 i 



4?r2 i(/i.2)+22^i(/i(«)/ 2 («),3) 
We evaluate the L-values that have popped up: 

= -£ x (z) - 7Bi(2«) + 8E 1 (4z) 
L(fx(u),s) = (-1 - 7 • 2-* + 8 • 4-*)C( S )L(x-3, «) 

L(/i(u),2) = (-1 - 1 + l)C(2)L(x-3,2) = -^ 2 L( X _ 3 ,2) 

and 

3 iq 

fi{u)fa{u) = --G 2 (z) - 5G 2 (2z) + yG 2 (3z) + 24G 2 (4z) - 35G 2 (6z) + 8G 2 (12z) 

3 19 

L(h(u)fa(u),s) = (- - - 5 • 2~ s + - ■ 3- s + 24 • 4" s - 35 • 6~ s + 8 • 12- s )C( S )C(s - 1) 

3 Fi 1Q 24 35 8 
L(/ 1 («)/ a (u),3) = (-3-^ + ^ + ^- ^+ iaS)C(3)C(2) 

= -H C(3 )C(2) = -^ 2 C(3) 
Eventually we finish the computation of the integrals in (|7.1[) : 

3 3 / 2 ; _ 7i 



^ 5i (*) • D- 2 g 2 (z)dz = -^L(x-3, 2) - — C(3), 



and multiplying by — 27rz we get 



^- 1 (. 9l ( <Z ) J D- 2 . 92 (g))|^ i = ^i( x _ 3 ,2) - |C(3) 



which is the same as (|6.4[) . 
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